Theory of Asymmetric Tunneling in the cuprate superconductors 
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We explain quantitatively, within the Gutzwiller-Resonating Valence Bond theory, the puzzling 
observation of tunneling conductivity between a metallic point and a cuprate high-T c superconductor 
which is markedly asymmetric between positive and negative voltage biases. The asymmetric part 
does not have a "coherence peak" but does show structure due to the gap. The fit to data is 
satisfactory within the over-simplifications of the theory; in particular, it explains the marked "peak- 
dip-hump" structure observed on the hole side and a number of other qualitative observations. This 
asymmetry is strong evidence for the projective nature of the ground state and hence for "t-J" 
physics. 
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In the conventional, BCS superconductors, the most 
complete and convincing evidence for the phonon mech- 
anism came from the tunneling spectrum-the tunnel- 
ing conductivity as a function of junction voltage. The 
features in this spectrum were shown to be uniquely 
caused by the anomalous self-energy (the "gap") and 
gave unequivocal evidence for its origin in the exchange 
of phonons. It has been a disappointment that tunneling 
spectroscopy has so far given us no such evidence for the 
cuprate superconductors. 

One of the puzzling experimental features about the 
tunneling in cuprates is the fact that the tunneling con- 
ductivity is markedly asymmetric as a function of volt- 
age. This was observed even in early crude attempts 
but became serious when vacuum tunneling from STM 
points achieved clean results in good detailed agreement 
(except for this fact) with the expected <i-wave density 
of states 0, . It is particularly interesting when one re- 
alizes that asymmetries are rare to non-existent in most 
metal-to-metal contacts, and are predicted not to exist 
(except for slow, continuous variations of tunneling prob- 
abilities) within Fermi liquid theory. The Gutzwiller- 
projected mean- field-theory is not a Fermi liquid 
theory and can - and does - exhibit asymmetry. 

The rarity of other examples of structure in tunneling- 
other than the well-known Giaevar effect of the appear- 
ance of the superconducting gap - is a consequence of 
two remarkable theorems. The first was proved by Har- 
rison [BJ essentially in order to explain why Giaevar did 
not see band structure effects in normal metals. Harri- 
son showed that the tunneling probability between two 
states, k in metal A and k' in metal B, evaluated in WKB 
approximation, contains a factor u(k)«(k') from the "at- 
tempt frequency" - where t>(k) is the velocity - in addi- 
tion to the WKB integral. (The theorem is actually more 
general than WKB, but this will do). This factor can- 
cels against the density of states, which is proportional to 
\/v. Although there may be prominent density-of-states 
anomalies in the band structure near to the Fermi level 
caused by a Van Hove singularity, they will not show up 
strongly in tunneling. 

The second theorem is due to Schrieffcr 6J. This is 
particularly useful in the BCS case, but has some bearing 
in the present one. Schrieffer pointed out that in most 



situations the tunneling probability is spread over a wide 
range of k-values, so that one must integrate the single- 
particle Green's function that appears in the tunneling 
conductivity over the variable k. This may be converted 
into a contour integral around the pole in the Green's 
function at the quasiparticle energy, and simplifies the 
tunneling density-of-states in the BCS case to 



N(E) = iV(0)Re • 



E 



^E 2 - A 2 (E) 



(1) 



Here A(E) is the gap function evaluated at the energy 
of the quasiparticle pole, where E 2 = [ek + S(k, E)] 2 + 
A(k. A';.-'. 

The result is not quite so clean in our case, where self- 
energies can be assumed to be k- as well as -E-dependent, 
but this modification seems only to require a factor of 
[1 + 9£/9ek] , which is likely to vary rather smoothly, 
at least near the coherence peaks. 

In the Gutzwiller mean-field-theory Q, we start from 
the approximation that the ground state of the t-J model 
Hamiltonian is a projected BCS product function, cho- 
sen by minimising the energy over all such functions. The 
t-J Hamiltonian is arrived at by a canonical transforma- 
tion of the true Hamiltonian, which presumably is 
essentially a Hubbard Hamiltonian: 

H tJ = e iS He- iS = PTP + J ^ S« ■ S 3 , (2) 

<i,j> 

where T is the kinetic energy and P the Gutzwiller pro- 
jection operator 



P = 



JJ(1 - n^mi) 



(3) 



The exchange term is not projected because it au- 
tomatically remains within the subspace defined by P, 
and hence commutes with it. The eigenstates of the t-J 
Hamiltonian are necessarily projected one-electron func- 
tions, so it is natural to approximate them with product 
functions. Thus the ground state variational function is 



I/) 



(4) 
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where 3>(A, fi) is the BCS product function 

m=Y[(u w + v k cl r cl kl \0), (5) 

k 

and Uk and are the variational parameters, determined 
by an effective BCS Hamiltonian that gives us a gap equa- 
tion as discussed in Ref. Q ■ This gap equation is actually 
the equation for the excitation energies of Gutzwillcr- 
projected excited-state wave functions 



|*k«r>=e< s £yLl*) 



where 



7 k | = u kC k | - •S'ukC-ki 



(6) 



(7) 



and the operator S creates a ground-state pair. The 
procedure is entirely analogous to the Hartree-Fock-BCS 
theory where the ground state is determined by the cri- 
terion that all single-Fermion excitations have positive 
energy. Within this theory, the excitations in Eq. H3 are 
the low energy single-Fermion excitations, by Koopman's 
theorem, which may be demonstrated in this case by the 
same method as in normal Hartree-Fock. 

We note that the theory so far, and its manipulations, 
are only correct because the Hamiltonian conserves parti- 
cle number, so that we do not need to consider coherence 
between states with different particle numbers. In order 
to specify that the number of electrons is correct we may 
simply fix average occupancies at the appropriate values, 
x for the empty state and — x) for the singly occupied 
ones of given spin; and then we proceed with Gutzwiller 
approximation based on those occupancies. But if, as 
in tunneling, we need to add or subtract electrons, we 
must follow Laughlin Q and introduce a fugacity factor 
Z for electron pairs. This is easily computed by noting 
that the ratio of these two occupancies in the original 
product function is (1 — x)/(l + x), so to get the cor- 
rect occupancies we must correct the normalization by 
(2x/(l + a;)) 1 / 2 . Hence for a pair of holes the fugacity 
factor is 



Z = 



2x 



(1 



(8) 



The resulting wave function, in the form given by Laugh- 
lin H, is Eq. multiplied by a factor £-(™T+ n i). (We 
choose for perspicuity to express Z as the fugacity of 
holes rather than electrons; the choice is arbitrary). 

This wave function may be rewritten in a form which 
demonstrates the effect of Z more clearly. In each factor 
(itk + WkC k jcL fc |) the i>k factor creates a pair of electrons, 
or conversely the u k factor can be taken as creating a 
pair of holes; thus in any component of the wave function 
which contains i*k, we can insert a corresponding factor 
Z. Thus instead of Eq. [5] we could write 



i$')=n(^+«k CkTC L ki )io), 



(9) 



and then we have absorbed the fugacity factor into $. 
But the individual factors in Eq. [5] are not normalized. 
To define the appropriate quasiparticle excitations as in 
Eq. they must be normalized and thus, finally, the 
appropriate product function must be written 



n 



(Zu k + ^kC kT c! ki ) 



u^Z 2 



10} 



(10) 



The individual factors define the product 7k7-k, so 
that the individual gamma contains the normalizing fac- 
tor (u 2 Z 2 + v 2 )* . Note that the inclusion of the Z factors 
nicely leads to the Z renormalization of the order param- 
eter, the superfluid density, and the kinetic energy 0. Z 
plays a role which can be described as the amplitude of 
the hole-pair wave function - or at least the relative am- 
plitude of the part of the pair function which is holes as 
opposed to spins. 

We would like to emphasize that the wave function (9) 
is completely identical to that used in previous papers 
and that no results as to quasiparticle energies or ground 
state averages are at all modified. The quantities u and 
v define the starting wave function which is to be pro- 
jected and its populations modified. We can think of the 
fugacity factor as part of the Gutzwiller projector, if we 
like; and like the projector, it does not commute with the 
fermion operators. Thus when we need to express the ex- 
citations in terms of real Fermions added to the system 
- outside the projector - we must use a modified set of 
u and v. But the real Fermions predominantly go in co- 
herently as single quasiparticles, except for the relatively 
small term for holes coming from the fluctuations of the 
opposite-spin occupancy. 

We want to present the simplest possible calculation, 
since the effect in question is a qualitative one. To this 
end we set e iS = 1, which involves errors of the order 
J/t which vary smoothly with energy - we are neglecting 
virtual double occupancy, and will therefore overestimate 
the asymmetry somewhat. 

Ignoring e , the quasiparticle wave function can be 
created from the unprojected product function |<f>") by 
either of the operators Pc^ or Pc_kj,. The former cre- 
ates it with relative amplitude uZj (Z 2 u 2 + v 2 ) i , the lat- 
ter with relative amplitude v/(Z 2 u 2 + v 2 )*. But what 
tunnels in from the metal is not a projected quasiparticle 
but a real one. We can think of the STM point as instan- 
taneously depositing a particle or hole into the Wannier 
function at the Cu atom under the point, and we then 
Fourier resolve the Wannier function amplitude at a time 
t = into excitations in the superconductor. Thus the 
operators which we must consider, acting on our assumed 
ground state, are cjP and CiP. For present purposes we 
can discuss only the site operators, later Fourier trans- 
forming to get the momentum space ones. 

Consider c\ a P. This may be divided into its two parts 
belonging to the forbidden and allowed subspaces: 



r P=(l-P)clP + PclP. 



(11) 
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The first term contains only excitations with energy 
larger than the Hubbard U and does not concern us. The 
second term is finite only x of the time; it requires that 
the site i contains an electron in the ground state function 
P|$"). Thus the probability that it creates an excitation 
in the allowed manifold is multiplied by x. But it will be 
important to note that because Pc'P = Per , obviously, 
the part of the hole that doesn't go into the upper Hub- 
bard band creates only single-quasiparticle excitations in 
this approximation, thus has only a coherent spectrum. 
Now consider cP. This automatically goes into the al- 
lowed manifold, but cP|$") is not exactly the same as 
the quasiparticle function Pc|<&") because the latter can 
contain components where |$") has m — 2 with proba- 
bility (1 — £ 2 )/4, while these do not appear in cP|$"). 
To adjust the normalization, note that 



c T P = c T (l - n T n x ) = P(l - n;)c T , 



(12) 



where the site index has been dropped. The average fac- 
tor reducing the quasiparticle function is 



cP = Pc(l- < n l >) 



(1 + x) 



u (hole-like) =v (electron- like), and vice versa. Thus the 
conductance is symmetric in the exchange of u and v, but 
NOT in the exchange of holes and electrons, in contrast 
to the BCS case. 

For electrons, the current is first of all reduced by the 
projection factor Z relative to that for the holes. Then 
the amplitude for a given channel is just the effective 
u for that channel, and the current its square; we get, 
taking into account the renormalization factor, that the 
tunneling density of states for electrons is 



Ne(E,A) = ^Z 



\Ju 2 + v 2 Z 2 \J v 2 + u 2 Z 2 



Pc. 



(13) 



(14) 

where u 2 = |[1 + e/E] and v 2 = |[1 - e/E]. 

For holes, we have no projection factor Z, but the hole 
amplitude contains the factor Z which can be thought of 
as the magnitude of the pair wave function. This satisfies 
the physical requirement that the coherent amplitude for 
holes and electrons must be the same at least at the same 
energy, and is also necessary for equilibrium. But the 
renormalization factor is not identical except at e = 0, 
E = A and rises as E — > oo to \/Z: 



Thus the ratio of the normalization factors for the elec- 
tron vs the hole spectra is gt — Z — 2x/(l + x). But in 
this case there is an incoherent spectrum, caused by the 
three- Fermion operator [n^— < n[ >]&; we do not be- 
lieve this is a large effect, but it may cause features in 
the hole spectrum, particularly in the neighborhood of A 
added to the magnetic resonance energy. 

For each k and spin, except for the rather small term 
mentioned in the last paragraph, there is only one quasi- 
particle wave function that appears in this approxima- 
tion, that obtained by Gutzwiller projecting the suitable 
BCS product function. Most of the amplitude is "coher- 
ent" , a conclusion quite different from that of Wen . 

At first sight one would think that the ratio of the 
tunneling conductivity for the sign of voltage V such as 
to inject a hole - external electrode positive - to that 
with the opposite sign would be just 2x/(l + x). But 
actually, quasiparticles are not pure electrons or holes 
but mixtures of the two, and precisely at the gap energy 
they are equal mixtures, so that at the gap energy the 
tunneling conductivity for +V and — V should be equal. 
The relevant tunnel current can flow either in the form 
of right-moving holes or left-moving electrons, and in the 
superconductor it is an equal coherent mixture of the 
two. But it is important to realize that for a given sign of 
voltage the two states which are coherent have actually 
the same charge, so that the hole is accompanied by a 
ground-state pair. 

Our calculation follows the method of Tinkham [Icj . 
As he points out (following in this Cohen, Falicov and 
Phillips H3) there are two quasiparticle channels with 
the same energy, with (e^ — Ef) positive and negative - 
electron-like and hole-like respectively. The itk and 
of the hole-like channel exchange their values, so that 



N h (E,A) = ^Z 



V u 2 + v 2 Z 2 \J v 2 + u 2 Z 2 



(15) 

Equations ^] and ^] show that the coherence peaks at 
e = are identical as predicted, but the ratio of the 
E — > oo asymptotes is Z , as expected from simple con- 
siderations. 
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FIG. 1: Tunneling conductance vs voltage for an optimally 
doped sample of BSCCO [data from S H Pan (unpublished)]. 
The cleanest data has been selected from among several sam- 
ples. Observations of J C Davis' group, and earlier work of 
Fischer 13), are similar in their general course. 

Finally it is necessary to take into account that we have 
a d-wave gap which means that the tunnel current must 
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be integrated over the gap distribution 
V{A)dA = "' A 



(16) 



with Ao the gap amplitude. 

The differential conductance for injection of electrons 
or holes is then 



G(E) t 



dA N e , h (E,A)V(A)dA, 



(17) 



with N et h given by Eqs. El an d El respectively. For 
numerical convergence, we add a weak imaginary part to 
the gap to simulate broadening, i.e. A — > (1 + if]) A with 
1) = 2 x 1CT 3 . 

We have approximated the Fermi Surface by a circle 
and normalized the maximum gap to 1. The result for 
the symmetric conductivity in the BCS case has often 
been displayed, and involves an elliptic function in its 
analytic form; there is a logarithmic peak at Ao, and a 
linear slope at E near zero. These shapes will appear 
in the region of the coherence peak and below in the 
Gutzwiller case, since the asymmetry in Eqs. 1141 and !T31 



is of greater than linear order in u 



e/E. But the 



renormalizations become appreciable even near E = Ao, 
and especially on the hole side can rise to dominate the 
peak for the underdoped case. The asymmetry has a 
pronounced upward cusp at E = Aq. It has been hard 
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FIG. 2: Computed curve of conductance vs voltage for Z 
= 0.2, which is too underdoped to fit Fig. (in which we 
estimate Z ~ 0.3) but the gap distribution is realistic to give 
a reasonable amplitude of coherence peak. 

to find an accurate curve for comparison of the predicted 
asymmetry with the published data. The best published 
data we can estimate from, for instance, Pan et al. 0, 
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FIG. 3: Computed curves for a sequence of Z's to demon- 
strate the variation with doping and to allow qualitative ex- 
trapolation to Z = 0. To achieve convergence, a complex 
gap A(l + ir/) (77 = 2 x 10 -3 ) is used in the distribution Eq. 
1161 The curves bear a rough resemblance to extremely under- 
doped experimental results. 



seems in surprisingly good agreement, especially when we 
realize that it is likely to contain traces of the interaction 
with the magnetic resonance of Keimer ^1 

The best curve from Pan's data is shown in Fig ^ 
which refers to a P^S^CaC^Os (BSCCO) sample at 
optimal doping. It should be noted that it is extremely 
difficult experimentally to make the conditions exactly 
identical for + and - V, and it is therefore likely that the 
very small asymmetry in the coherence peak structure is 
an experimental artifact, and is actually absent. Ignor- 
ing that, the fit to the general course of the asymmetry 
is remarkable. Other data not shown here confirm the 
general behavior with doping as well. 

A better estimate of the experimental curve would take 
into account that the band is much flatter at the gap 
antinodes near 0, so that the maximum of the gap should 
be more strongly emphasized. An estimate of such an ef- 
fect can be obtained by adding to V(A) in Eq. El a 
delta- function of unity at A ; Fig. [3 shows the result of 
such an exercise, which seems to reproduce the experi- 
mental curve of Pan at optimal doping fairly well. 

Another somewhat speculative exercise is to continue 
Z towards 0, which gives us a conjectured tunneling spec- 
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trum for the pure RVB phase, which is our model for the 
pseudogap state, at least at higher temperatures. In Fig. 
we show how the variation with doping goes. In the 
limit as Z — > 0, the current (almost all on the hole side) 
does not extrapolate to an asymptote but continues to 
rise linearly at high voltage. The ratio of the asymptotes 
on the two sides is 1/Z, an observation which seems to 
accord with most estimates of doping percentages. 

The curves which represent regions the experimental- 
ists think are quite underdoped differ from higher dopings 
in that the symmetrical parts of the curve extend only 
to rather low voltages, and the coherence peaks are sup- 
pressed. The higher-voltage conductivity seems almost 
completely composed of the asymmetric "hump" behav- 



ior and to be dominantly on the hole side. In the same 
regions, a characteristic "4 x 4" density wave is observed. 
In a forthcoming paper we will suggest a mechanism that 
might relatively suppress the gap antinodes. 
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